Interacting spinor and scalar fields in Bianchi type-I Universe filled with 
viscous fluid: exact and numerical solutions 

Bijan Saha 

Laboratory of Information Technologies 
Joint Institute for Nuclear Research, Dubna 
141980 Dubna, Moscow region, Russia^ 
(Dated: September 7, 2008) 

We consider a self-consistent system of spinor and scalar fields within the framework of a 
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I. INTRODUCTION 

The investigation of relativistic cosmological models usually has the energy momentum tensor 
of matter generated by a perfect fluid. To consider more realistic models one must take into 
account the viscosity mechanisms, which have already attracted the attention of many researchers 

aaaaaaaaaGa. 

The nature of cosmological solutions for homogeneous Bianchi type I (BI) model was investi- 



gated by Belinsky and Khalatnikov ||11|1 by taking into account dissipative process due to viscosity. 
In [fl^.ll3ll we reinvestigate the problem posed in [[Till in presence of a A term. Though that Mur- 
phy [9] claimed that the introduction of bulk viscosity can avoid the initial singularity at finite 
past, but Belinsky and Khalatnikov 111 11 1 showed that viscosity cannot remove the cosmological 
singularity but results in a qualitatively new behavior of the solutions near singularity. To elimi- 
nate the initial singularities a self-consistent system of nonlinear spinor and BI gravitational field 
was considered by us in a series of papers [Hi [3 Oil For some cases we were able to find 



field configurations those were always regular. Recently it was found that the introduction of a 
spinor field into the system may explain the late time acceleration of the Universe, hence can be 
considered as an alternative to dark energy fl 1 3M - 

Given the importance of viscous fluid and spinor field to construct a more realistic model of the 
Universe, we in ||19Ll20|l we introduced spinor field into the system and solved the system for some 
special choice of viscosity. The purpose of this paper is to study an interacting system of spinor 
and scalar fields within the scope of a Bianchi type I cosmological model filled with viscous fluid 
in presence of a A term and clarify the role of viscosity and field interaction in the evolution of the 
Universe. 



*Electtonic address: saha@thsunl.jinr.ru, bijan@jinr.ru UPvL: http://thsunl.jinr.ru/~sciha/ 
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II. DERIVATION OF BASIC EQUATIONS 

In this section we derive the fundamental equations for the interacting spinor, scalar and gravi- 
tational fields from the action and write their solutions in term of the volume scale T defined bellow 
(12.161) . We also derive the equation for % which plays the central role here. 

We consider a system of nonlinear spinor, scalar and BI gravitational field in presence of perfect 
fluid given by the action 

^(g;y,xjr) = J J?^gd£l (2.1) 

with 

= =S^r 4- =S?ss ~t~ =^^m- (2.2) 

The gravitational part of the Lagrangian (12.21) is given by a Bianchi type I (BI hereafter) space- 
time, whereas «if ss describes the interacting spinor and scalar field lagrangian and Jz? m stands for 
the lagrangian density of viscous fluid. 



A. Material field Lagrangian 



We choose the interacting spinor and scalar field Lagrangian as 



(2.3) 



Here m is the spinor mass, X is the coupling constant and F = F(I,J) with I = 1$ = S 2 = (yy) 2 
and J = Ip = P 2 = (zty/y 5 ty/) 2 . We would like to mention that there are 5 invariants constructed from 
bilinear spinor from. Using Fierz transformation it can be shown that among the five invariants 
only I and J are independent as all other can be expressed by them: Iy = —Ia = h + h and Iq — 
2 (7s — Ip) 112 U l22l |23L |24|]. Since the the bilinear identities appear to have been given in literature 
only partially and the relation between the invariants are given as problem [cf. eg. [|24l |25fl ], we 
work them out in the appendix below. Therefore, the choice F = F(I,J), describes the nonlinearity 
in the most general of its form 111511 . Note that setting A = in (|2.3I) we come to the case with 
minimal coupling. 



B. The gravitational field 

As a gravitational field we consider the Bianchi type I (BI) cosmological model. It is the 
simplest model of anisotropic universe that describes a homogeneous and spatially flat space-time 
and if filled with perfect fluid with the equation of state p = £e, £ < 1, it eventually evolves 



into a FRW universe H261 . 12711 . The isotropy of present-day universe makes BI model a prime 
candidate for studying the possible effects of an anisotropy in the early universe on modern-day 
data observations. In view of what has been mentioned above we choose the gravitational part of 
the Lagrangian (|2.2I) in the form 

^ B = 4 (2-4) 

where R is the scalar curvature, K = SnG being the Einstein's gravitational constant. The gravita- 
tional field in our case is given by a Bianchi type I (BI) metric 

ds 2 = dt 2 -a 2 dx 2 -b 2 dy 2 -c 2 dz 2 , (2.5) 



with a, b, c being the functions of time t only. Here the speed of light is taken to be unity. 
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C. Field equations 

Let us now write the field equations corresponding to the action (12.11) . 
Variation of (12.11) with respect to spinor field yf(\fr) gives spinor field equations 

*'7 M V M v/-mv/+^v/ + ^?'yV = 0, (2.6a) 

Np\j/yV+m\j/-@\j/-&i\j/y 5 = 0, (2.6b) 

where we denote 

_ A a dF X a dF 

Varying (|2.1I) with respect to scalar field we find 

1 d 



(v / =gg VAi (l+^)(p, M )=0. (2.7) 





(2.8a) 


kT^ + A, 


(2.8b) 


kT^+A, 


(2.8c) 




(2.8d) 



V=gdx v 

Variation of (|2.1I) with respect to metric tensor g^ v gives the Einstein's field equation. For the 
BI space-time (12.51 ) on account of the A term this system has the form 

b c be 

T+-+T- 

b c be 

c a c d 

- + -H 

c a c a 

a b db 

a b ab 
db be cd 
ab be ca 

where over dot means differentiation with respect to t and Ty is the energy-momentum tensor of 
the material field given by 

t£ = ^{wy^vW+wyv^iiW-^^wyvw-^vwy^ (2.9) 

+ (\-XF) W P-8Pj? + T jl v m . 
Here T^ m is the energy-momentum tensor of a viscous fluid having the form 

TjTm = ( e + p')u^u v - + r]g yl3 [m m;J3 + up.^ - u^u a up. a - upu a u^ a ] , (2.10) 

where 

p'=P-^-\n)^- (2.1D 

Here e is the energy density, p - pressure, r\ and ^ are the coefficients of shear and bulk viscosity, 
respectively. In a comoving system of reference such that = (1, 0, 0, 0) we have 

r ° m = £, (2.12a) 
Tim = ~P' + 2 ^ < 2 - 12b ) 

TL = -p' + 2r]~, (2.12c) 
Tj m = -p' + 2r] C -. (2.12d) 
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In the Eqs. (12.61) and (12.91 ) is the covariant derivatives acting on a spinor field as 11281. 12911 



dw d\if 

v MV = ^r-r M v, v M ?=^+yr M> (2.13) 



where are the Fock-Ivanenko spinor connection coefficients defined by 

r M = Jy ff (rJ arv -a M y CT ). (2.14) 
For the metric (12.51) one has the following components of the spinor connection coefficients 

r = 0, T x = l ~a{t)ff, T 2 = ^b(t)ff 1 T 3 = l -c{t)ff. (2.15) 

The Dirac matrices y" (x) of curved space-time are connected with those of Minkowski one as 
follows: 

7° = 7°, Y l = ?/a, f = f/b, y i = f/c 



with 



^<o-°,y ?=(-A\ 



where <x are the Pauli matrices: 



„> = r?n. o*= °.-J . o 3 " 



i o y ' y i / ' v o — l 

Note that the f and the o matrices obey the following properties: 

ffJ + fJf = 2r 1 i ^ ij = 0,1,2,3 
ff + ff = 0, [ff=I, i = 0,l,2,3 
ct'ct* = + /e jW a z , j,k,l = 1,2,3 

where T] ;/ = {1,-1,-1,-1} is the diagonal matrix, is the Kronekar symbol and e^jy is the 
totally antisymmetric matrix with £123 = +1. 

We study the space-independent solutions to the spinor and scalar field equations (12.61) so that 
y = yr(t) and <p = <p(t). Here we define 

x = abc=y^g (2.16) 
Under this assumption from (|2.7I) for the scalr field we find 

<p = C J[x{l+XF)Y l dt, C = const. (2.17) 

The spinor field equation (I2.6al) in account of (|2.13l) and (12.151) takes the form 

/y°( 4- + -^ ) V/-mv/ + ^v/ + ^z'y 5 v/ = 0. (2.18) 
\ot 2x J 
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Setting Vj(t) = y/xy/jit), j = 1,2,3,4, from (12.181) one deduces the following system of equa- 
tions: 

Vi+i(m-Si)Vi-9Vs = 0, (2.19a) 

V 2 + i{m-S>)V 2 -^V4 = 0, (2.19b) 

%-i{m-@)V$+9V\ = 0, (2.19c) 

V 4 -i(m-@)V 4 + yV2 = 0. (2.19d) 

From (|2.6al ) we also write the equations for the invariants 5, P and A = y/ v j/ 5 y°y/' 

S -2^A Q = 0, (2.20a) 

P -2(m-@)A = 0, (2.20b) 

A + 2(m-3!)Po + 2&So = 0, (2.20c) 

where So = tS, Po = TP, and Aq = xA. The Eq. (|2.20l) leads to the following relation 

S 2 + P 2 +A 2 = cl/z 2 , C\ = const. (2.21) 

Giving the concrete form of F from (12.191 ) one writes the components of the spinor functions 
in explicitly and using the solutions obtained one can write the components of spinor current: 

f = ijrfy. (2.22) 

The component j° 

f = - \v{Vi + V 2 *V 2 + V 3 * V 3 + V 4 *V 4 ] , (2.23) 

T 

defines the charge density of spinor field that has the following chronometric-invariant form 

P = 0W ) 1/2 (2-24) 
The total charge of spinor field is defined as 

oo 

Q= J p^- r gdxdydz = pxt, (2.25) 

— oo 

where Y is the volume. From the spin tensor 

S HV,e = 1 ytfo^ + O^fty. (2.26) 

one finds chronometric invariant spin tensor 

4f=(%0^°) 1/2 , (2.27) 
and the projection of the spin vector on k axis 

s k= J S i ^^/- r gdxdydz = S i ^xV. (2.28) 
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Let us now solve the Einstein equations. To do it we first write the expressions for the compo- 
nents of the energy-momentum tensor explicitly: 



T n u = mS + 



C 2 



2x 2 (l+XF) 
C 2 



+ £ 



2t 2 (1+AF) 



&S + &P 



&S + &P 



C 2 



-p' + 2T}- 
a 

b 



a 



2x 2 (l+XF) 
C 2 

2t 2 (1+AF) 



- p + 2t) 



P+2t) 



Tl+27] 



(2.29a) 
(2.29b) 
(2.29c) 
(2.29d) 



In account of (12.291) from (12.81) we find the metric functions [15] 

e -2KjT]dt 



a(t) 
b{t) 
c(t) 



yiT 1/3 exp 
y 2 T 1/3 exp 
F 3 T 1/3 exp 



3 
3 

X3 
3 



T(0 
e -2x fj}dt 

e -2K jr\dt 



no 



(2.30a) 
(2.30b) 
(2.30c) 



with the constants Y[ and X[ obeying 

Y 1 Y 2 Y 3 = l, 



X 1 +X 2 + X 3 = 0. 



As one sees from (|2.30al) . (|2.30bl) and (|2.30cl) . for x = t" with n > 1 the exponent tends to unity 
at large t, and the anisotropic model becomes isotropic one. 

Further we will investigate the existence of singularity (singular point) of the gravitational 
case, which can be done by investigating the invariant characteristics of the space-time. In general 
relativity these invariants are composed from the curvature tensor and the metric one. In a 4D 
Riemann space-time there are 14 independent invariants. Instead of analyzing all 14 invariants, 
one can confine this study only in 3, namely the scalar curvature I\ = R, I% = R^yjlV, and the 

Kretschmann scalar I 3 = R a p^iv^ a ^ V ■ At any regular space-time point, these three invariants 
^l; h, h should be finite. One can easily verify that 



h 



h 



h 



Thus we see that at any space-time point, where x = the invariants I\, I2, h, as well as the scalar 
and spinor fields become infinity, hence the space-time becomes singular at this point. 

In what follows, we write the equation for x and study it in details. 

Summation of Einstein equations (I2.8al) . (I2.8bl) . (I2.8cl) and (I2.8dl) multiplied by 3 gives 



f = ~K(T£ + T l l )x + 3K7ix + 3Ax, 



which can be rearranged as 



t-^k\!;t = ^ic(mS + @S + &P + £ - p^Jx + 3Ax. 



(2.31) 



(2.32) 
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For the right-hand-side of (|2.32l) to be a function of x only, the solution to this equation is well- 
known B0I1 . 

On the other hand from Bianchi identity G^. v = one finds 

7m v ; v = Tv v +r v pv T£ -rj v r p v = o, (2.33) 

which in our case has the form 

I( Tr oy _ * T } _ K2 _ c T 3 = Q (2 34) 

*V (X o c 

This equation can be rewritten as 

Recall that (12.201) gives 

(m-@)S -&P = 0. 
In view of that after a little manipulation from (12.351) we obtain 

£ + -co - (£ + Itj ) ^ + 4?7 (kTo + A) = 0, (2.36) 

T 3 T z 

where 

u) = £ + p, (2.37) 

is the thermal function. Let us now in analogy with Hubble constant introduce the quantity H, 
such that 

x d b c 

- = - + - + - = 3if. (2.38) 
X a b c 

Then (12.321) and (12.361) in account of (12.291) can be rewritten as 

H = | (3%H - a) - (3H 2 - ke - A) + | (mS + ®S + 9P) , (2.39a) 



£ = 3H(3%H-(d) +4ri(3H 2 -K£-A) -ArjK 



mS + 



C 1 



2t 2 (1+AF)J 



(2.39b) 



Thus, the metric functions are found explicitly in terms of X and viscosity. To write x and compo- 
nents of spinor field as well and scalar one we have to specify the function F. In the next section 
we explicitly solve Eqs. (12.191) and (12.391) for some concrete value of F. 



III. SOME SPECIAL SOLUTIONS 

In this section we first solve the spinor field equations for some special choice of F, which will 
be given in terms of x. Thereafter, we will study the system (12.391 ) in details and give explicit 
solution for some special cases. 



A. Solutions to the spinor field equations 



As one sees, introduction of viscous fluid has no direct effect on the system of spinor field 
equations (12.191) . Viscous fluid has an implicit influence on the system through x. A detailed 
analysis of the system in question can be found in IU5I1 . Here we just write the final results. 
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1. Case with F = F{I) 

Here we consider the case when the nonlinear spinor field is given by F = F(I). As in the case 
with minimal coupling from (12.20al) one finds 

S=— , Co = const. (3.1) 

T 



For components of spinor field we find 111511 



V T V T 

(3.2) 

V T 

with C, being the integration constants and are related to Co as Co = C\ + C\ — C| — C|. Here 
P = f(m-@)dt. 

For the components of the spin current from (12.221) we find 

f = ~[C? + c£ + Cf + c2], j 1 = ^[C 1 C 4 + C 2 C 3 ]cos(2/3), 
/ = A[ Cl C 4 -C 2 C 3 ] sin(2j8), j 3 = -[C 1 C 3 -C 2 C 4 ]cos(2 J 8), 
whereas, for the projection of spin vectors on the X, Y and Z axis we find 

^,23,0 _ C1C2+C3C4 ^ 31)0 _ q ^,12,0 _ C l ~ ^2 + ^3 ~ ^4 

Total charge of the system in a volume ^ in this case is 

Q = [C? + Cf + C 3 2 + Cf ] r. (3.3) 

Thus, for t 7^ the components of spin current and the projection of spin vectors are singularity- 
free and the total charge of the system in a finite volume is always finite. Note that, setting A = 0, 
i.e., j8 = mt in the foregoing expressions one get the results for the linear spinor field. 

2. Case with F = F {J) 

Here we consider the case withF =F(J). In this case we assume the spinor field to be massless. 
Note that, in the unified nonlinear spinor theory of Heisenberg, the massive term remains absent, 
and according to Heisenberg, the particle mass should be obtained as a result of quantization of 



spinor prematter p In the nonlinear generalization of classical field equations, the massive 
term does not possess the significance that it possesses in the linear one, as it by no means defines 
total energy (or mass) of the nonlinear field system. Thus without losing the generality we can 
consider massless spinor field putting m = 0. Then from (12.20bl) one gets 



P = Dq/%, Dq = const. 



(3.4) 
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In this case the spinor field components take the form 

V X v T 

(3.5) 

The integration constants D, are connected to Do by Do = 2 (D 2 + F>\ — D 2 — D 2 ) . Here we set 
o = f&dt. 

For the components of the spin current from (12.221) we find 

f = — + D 2 + D 3 + D4] , / = — [D 2 D 3 +D 1 D 4 ]cos(2a), 
j 2 = y % [D2D3 -D1D4] sin(2a), / = — [D1D3 — D2D4] cos(2ct), 
whereas, for the projection of spin vectors on the X, Y and Z axis we find 

523i o_ 2(D 1 D 2 +D 3 D 4 ) ^ S 3l,0 = 0j s i2.o_ Dj-D 2 2+D 2 3-Dl 
hex 2abx 

We see that for any nontrivial x as in previous case the components of spin current and the pro- 
jection of spin vectors are singularity-free and the total charge of the system in a finite volume is 
always finite. 



B. Determination of x 

In this subsection we simultaneously solve the system of equations for x and e. Since setting 
m = in the equations for F = F(I) one comes to the case when F = F(J), we consider the case 
with F being the function of / only. Let F be the power function of 5, i.e., F = S". As it was 
established earlier, in this case 5 = Cq/x, or setting Co = 1 simply S = l/x. for simplicity we also 
set C = 1 . Evaluating @ in terms of X we then come to the following system of equations 



r n-l 



3 „ . 3k Xn x' 

r + TlT + T(A + t«) 2 



m 



L X 



+ 



3k ( 

% n-2 

2(X + x n ) 



p)x- 3At, 



+ A 



or in terms of H 



(3.6a) 
(3.6b) 



x 
H 



3Hx 
1 



(3^H-0))-(3H 2 -e-A) + 



K \m Xn x n 2 
X + T(A+T' 3 ) 2 



e = 3H(3$H- go) +4r}(3H 2 -£-A) -Ar\ 



m x n 2 
T + 2(A + t") 



(3.7a) 
(3.7b) 

(3.7c) 



Here 7] and are the bulk and shear viscosity, respectively and they are both positively definite, 
i.e., 

r\ > 0, % > 0. (3.8) 
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They may be either constant or function of time or energy. We consider the case when 

T] =Ae a , ^=Be p , (3.9) 

with A and B being some positive quantities. For p we set as in perfect fluid, 

P = Ce, Ce(0,l]. (3.10) 

Note that in this case £ ^ 0, since for dust pressure, hence temperature is zero, that results in 
vanishing viscosity. 

The system (13.71 ) without spinor field have been extensively studied in literature either partially 
or as a whole fl 1 IN . Here we try to solve the system (|3.6I) for some particular choice of 



parameters. 



1. Case with bulk viscosity 



Let us first consider the case with bulk viscosity alone setting coefficient of shear viscosity 
77 = 0. In this case from (12.311) and (12.351) we find the following relation 

jc? ° = 3H 2 - A + C 00 , Coo = const. (3.11) 

We also demand the coefficient of bulk viscosity be inverse proportional to expansion, i.e., 

^6 = 3^H = C 2 , C 2 = const. (3.12) 

Inserting r\ = 0, (T3.12I) and (13.101) into (|3.7cl) one finds 

£= T ^[C 2 -C 3 /x 1+ Z}. (3.13) 

Then from (|3.6a|) we get the following equation for determining t: 



3 



Km + 3 



2 L 2 



C2 K + A 



3k(1-C)C 2 t 1+ £-C3 3kXh x n ~ l . ,„ i/lx 

T+ 2(1 + Tt - + — (UW E%t) ' <3 ' 14) 



where g is the set of problem parameters. As one sees, the right hand side of the Eq. (13.141) is 
a function of x, hence can be solved in quadrature [|3Qh . We solve the Eq. (13.141) numerically. It 
can be noted that the Eq. (13.141) can be viewed as one describing the motion of a single particle. 
Sometimes it is useful to plot the potential of the corresponding equation which in this case is 

W(q,r) = -2 J &{q,x)dx. (3.15) 

The problem parameters are chosen as follows: K = 1, m = 1, A = 0.5, £ = 1/3, n = 4, C 2 = 2 and 
C3 = 1. Here we consider the cases with different A, namely with A = —2,0, 1, respectively. The 
initial value of x is taken to be a small one, whereas, the first derivative of x, i.e., x at that point 
of time is calculated from (13.111) . In Fig. [T]we have illustrated the potential corresponding to Eq. 
(13.141) . It can be immediately seen that independent of the sign of A we have always expanding 
universe. But as is seen from Fig. [2] a positive A results in accelerated mode of expansion, while 
the negative one causes deceleration. 
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FIG. 1 : View of the potential corresponding to 
the different sign of the A term. 



FIG. 2: Evolution of z depending on the signs of 
the A term. 



2. Case with bulk and shear viscosities 



Let us consider more general case. Following [12] we choose the shear viscosity being propor- 
tional to the expansion, namely, 

3 1 

ri = -—H = -—e. (3.16) 



2k 2k 
In absence of spinor field this assumption leads to 



3H 2 



KE+C4, C4 = const. 



(3.17) 



It can be shown that the relation (13.171) in our case can be achieved only for massless spinor field 
with the nonlinear term being 

F=(F S 2 -1)/X. 



Equation for % in this case has the form 



Tf -0.5(1 - Qi 2 - 1.5k%tt -3[A-0.5(1 - QC 4 + k/2F ]t 2 = 0. 



(3.18) 



In case of t, = const, there exists several special solutions available in handbooks on differential 
equations. For that reason we rewrite this equation in terms of H: 



H = -1.5(1 + QH 2 + \.5k$H+ [A- 0.5(1 - QC 4 + k/2F }. 
The solution of the foregoing equation can be written in quadrature as 

dH 

AH 2 + BH + C 



(3.19) 



(3.20) 
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withA = -1.5(l + £),fl = 1.5k<§ andC = A-0.5(l-C)C 4 + fc/2F . If the bulk viscosity is taken 
to be a constant one, i.e., £ = const., then depending on the value of the discriminant B 2 — AAC 
there exists three types of solutions, namely 113411 : 



^y Mg+i+^ i B 2 >AAC, 



VS 2 -4AC 1 2A//+B-Vfi 2 -4AC 



V4AC-S 2 



arctan 



2AH+B 
V4AC-B 2 ' 



2AH+B1 



B 2 < AAC, 



B 2 = AAC. 



(3.21) 



Note that a detailed analysis of these solutions in absence of spinor and scalar field was given 
in fl2|l . We choose the problem parameters as follows: C4 = 9, £ = 1/3 and K" = 2. Under this 
choice we find 8 = B 2 — AAC = 9^ 2 + 8 ( A — 2) . After that we chose A positive, trivial or negative 
(in particular we chose A = (7/8,0, —5/2)). The quantity £ now is taken such a way that we have 
S = (Si, &i, £3) for all values of A chosen above, whereas, Si > 0, &i — and £3 < 0. In Figs. 
[3l |4] and [5] we plot the evolution of T corresponding to a trivial, positive and negative value of S, 
respectively. As is seen, the behavior of X mainly depends on S and independent to the sign of A. 
Since a negative 8 gives non-periodic mode of evolution we plot the corresponding phase diagram 
in Fig. [6] 





FIG. 3: Evolution of the universe with 8 = for 
different A. 



FIG. 4: Evolution of the universe with 8 > for 
different A. 



IV. CONCLUSION 

We consider the self consistent system of spinor, scalar and gravitational fields within the 
framework of Bianchi type-I cosmological model filled with viscous fluid. Solutions to the cor- 
responding equations are given in terms of the volume scale of the BI space-time, i.e., in terms 
of T = abc. The system of equations for determining T, energy-density of the viscous fluid £ and 
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A = 



A>0 





FIG. 5: Evolution of the universe with 8 < for 
different A. 



FIG. 6: Phase diagram corresponding to the case 
with 8 < for different A. 



Hubble parameter H has been worked out. Exact solution to the aforementioned system has been 
given only for some special choice of viscosity. It should be noted that the system (13.71) is far richer 
and allows a number of mathematically interesting results, though not all of them is physically re- 
alizable. Given this fact we plan to review this system and give a detailed analysis and qualitative 
solutions of the corresponding system in some of our future works. 

V. APPENDIX 



Let us now construct the invariants of spinor field. Since \\f and yr* (complex conjugate of 
has 4 component each, one can construct 4 • 4 = 16 independent bilinear combinations. They are 
scalar, pseudoscalar, vector, axial vector, and tensor denoted, respective, by 



S 
P 

2 MV 



(5.1a) 
(5.1b) 
(5.1c) 
(5.1d) 
(5.1e) 



where o^ v = (z'/2)[y^y v — y v y^] is the anti- symmetric tensor. Invariants, corresponding to the 
bilinear forms are 

S 2 = (W) 2 , (5.2a) 

P 2 = (zV/yV) 2 , (5.2b) 



Is 
Ip 
h 
h 



(5.2c) 
(5.2d) 
(5.2e) 
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/matrices in the above expressions obey the following algebra 

and are connected with the flat space-time Dirac matrices f in the following way 

8fiv{x) = 4W4W^> 7m W = 

where T]^ = diag( 1,-1,-1,-1) and is a set of tetrad 4- vectors. 
For the diagonal metric such as Bianchi-I 



we have 



ds 2 = dt 2 - a 2 (t)dx 2 - b 2 (t)dy 2 - c 2 (t)dz 2 



To = % Y\=a{t)% J2 = b{t)% y i = c {t)% 



7° = 7°, 7 1 = 7 1 M0, f = f/Ht), 7 i = f/c(t). 
Flat space-time matrices f we will choose in the form: 

/ 1 



7° 



f 



Defining y 5 as follows, 



/l \ 
10 0' 
0-10 

Vo o o -l / 

/ -i\ 
i 
0/00 

V -i o o o / 



7 



T 3 



10 
0-100 

V -l o oo 

/ 1 \ 
-1 
-10 

V o io o / 



T 5 = -tWtVtV, £p 



^-p.upi t i t i -^pivap — V~g £ y.v<yp, £0123 — 1? 



we obtain 



f 



( -1 \ 
-1 
-10 

V o -l o o / 



Note that y is a 4 component function given by, 



¥2 
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Denoting F bilinear spinor form in Minkowski spacetime and F in curve spacetime (in our case in 
BI) we find the following expressions for the non-trivial components of bilinear spinor form: 



S = ( y/JYi + x\r 2 - ¥3 ¥3 - Wt Va), S = S, (5.7) 
P = -i{Yi ¥3 + ¥2 ¥4 - ¥3 ¥\ ~ ¥a ¥2), P = P, (5.8) 
^ = (^1 + ^ + ^3 + ^4), V° = V°, (5.9) 

V 1 = (^i>4+^2>3 + ^3>2+^4>l), V l =V l /a, (5.10) 

^ = -1(^4-^3 + ^-^1), V 2 = V 2 /b, (5.11) 

^ = (^3-^4 + ^1-^), V 3 =V 3 /c, (5.12) 

A = (¥t¥3 + ¥^¥4 + ¥l¥i + ¥4¥2), A°=A°, (5.13) 

A 1 = (¥t¥2 + ¥2*¥i + ¥3¥4 + ¥4¥3), A l =A l /a : (5.14) 

A 2 = -/(WVi-V^Vi + V^V^-V^V^). A 2 =A 2 /b, (5.15) 

A 3 = (^i>i-V^2>2 + ^3>3- V4>4), A 3 =A 3 /c, (5.16) 

<2 01 = /(^> 4 + ^2>3-^3>2-^4>l), Q 0l =Q° l /a, (5.17) 

2 02 = (^> 4 -^2>3-^3>2+^4>l), Q 02 = Q° 2 /b, (5.18) 

Q 03 = /(^>3-^2>4-%>l + ^4>2), Q 03 = Q° 3 /C, (5.19) 

<2 12 = (^i>i-^2>2-^3>3 + ^4>4), Q l2 = Q l2 /ab, (5.20) 

2 23 = (^> 2 + ^2>l-^3>4-^4>3), Q 23 = Q 23 /bc, (5.21) 

e 13 = /(^> 2 - v^2>i - %>4 + y 4 >3), e 13 = e 13 /«c. (5.22) 

(5.23) 

Using the above expressions we find 

I S = 5 2 = (v/ 1 V 1 ) 2 +(v/ 2 > 2 ) 2 +(^3V3) 2 + (^4>4) 2 

+ 2[Yt ¥1 ¥2 ¥2 - ¥1 ¥1 ¥3 ¥3 ~ ¥1 ¥1 ¥l ¥4 

- ¥l¥2¥l¥3-¥l¥2¥4*¥4 + ¥l¥3¥l¥4]- (5.24) 

I P = J P 2 = -( V , 1 V 3 ) 2 -( V /| V / 4 ) 2 -( V ,| V / 1 ) 2 -(^4V2) 2 

- 2[i^i ¥3 ¥2 ¥4 - ¥1 ¥1 ¥3 ¥3 ~ ¥1 ¥3 ¥4 ¥2 

- ¥2*¥4¥l¥\-¥2*¥2¥Z¥4 + ¥3¥i¥l¥2]- (5.25) 



Iy = (V ) 2 -(V l ) 2 -(V 2 ) 2 -(V 3 ) 2 

= (¥1 ¥i) 2 + (¥2 ¥2? + (¥3 ¥3? + (¥l ¥4) 2 

- (¥t¥3) 2 -( ¥l ¥4) 2 -(¥l¥i) 2 -( ¥l ¥2? 

+ 2 [ y{ xj/x y/2 - V^i V 1 ¥4 ¥4 - ¥2 ¥2 ¥3 ¥3 + ¥3 ¥3 ¥4 ¥4 

- ¥1 ¥4 ¥2 ¥3 ~ ¥3 ¥2 ¥l ¥1 + ¥1 ¥3 ¥4 ¥2 + ¥2 ¥4 ¥3 ¥1 ] 

= Is + Ip- (5.26) 
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I A = (A ) 2 — (A 1 ) 2 — (A 2 ) 2 — (A 3 ) 2 

= -(Vi>i) 2 - (V2 V2) 2 - (¥l¥3) 2 - (¥l¥4) 2 

+ ( W Vs) 2 + ( Vl Wa) 2 + ( Wl¥i f + ( ¥4 W2) 2 

- 2[\j/l y/i y/| y/2 - K y/4* y/ 4 - y/f y^y! ¥3 + ¥3 ¥3 ¥l ¥4 

- ¥1 ¥4 ¥2 ¥3 - ¥3 ¥2 ¥4 ¥1 + ¥1 ¥3 ¥l ¥2 + ¥2 ¥4 ¥3 ¥1 } 

= -(Is + Ip). (5.27) 



Iq = 2[(Q 12 ) 2 + (Q 13 ) 2 + (Q 23 ) 2 - (Q 01 ) 2 - (Q 02 ) 2 - (Q 03 ) 2 ] 
= (¥l¥i) 2 + (¥2¥2) 2 + (¥3¥3) 2 + (¥4 ¥ ¥4) 2 
+ (¥t¥3) 2 + (¥2 f ¥4) 2 + (¥l¥i) 2 + (¥l¥2) 2 
+ 2[y4 Yi ¥2 ¥2 - ¥1 ¥1 ¥4 ¥4 - ¥2 ¥2 ¥3 ¥3 + ¥3 ¥3 ¥4 ¥4 
+ \i4 y/ 4 y/| y/ 3 + ¥3 ¥2 ¥4 ¥1 ~ ¥1 ¥3 ¥4 ¥2 - ¥2 ¥4 ¥3 ¥1 } 
- 4(v/ 1 > 1 v/ 3 >3 + Y/|y/2V 4 >4) = 2(Is-I F ). (5.28) 

Thus we see that that the invariants I v , I A and Iq can be expressed in terms of Is and Ip. 

An alternative proof of Iy = —I A = (Is +Ip) and Iq = 2(Is — Ip) can be given using Feirz trans- 
formation. In doing so we recall that any 4x4 matrix T can be presented as a linear combination of 
/: 

r = £c A /\ c A = -Try A Y. (5.29) 

A 4 

Explicitly it can be presented as 

r ik = W r imtiirAik- (5.30) 

4 A 

Here 

Y A = {I,Y 5 ,f,if"f,i°^ V }, A=l,2..16. (5.31) 
The completeness condition gives 

¥« = !l^'- (532) 

4 A 

Multiplying the foregoing expression by ¥i¥k¥m¥i one g ets 

(vVXvV) = j J L(rrA¥ b )(VY¥ cl )- (5.33) 

4 A 

Replacing y/^ — > y 8 ^ 1 and yr — > "f\\r and using the expression 

rV = LQ/, C R = 1 -Try A fy R , (5.34) 
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one gets the other identities . Denoting I s = ( y" y b ) ( \f ^ ) , l' s = ( ^ ) ( \f y b )•••■■ ■ one finds 



After a little manipulations from the foregoing system one comes to desired relations between the 
invariants of the bilinear spinor fields. 
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